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We report the electric transport study of the three-dimensional topological insulator TlBiSe 2 • We 
applied a newly developed analysis procedure and precisely determined two-carrier transport prop¬ 
erties. Magnetotransport properties revealed a multicarrier conduction of high- and low-mobility 
electrons in the bulk, which was in qualitative agreement with angle-resolved photoemission re¬ 
sults [K. Kuroda et al., Phys. Rev. Lett. 105, 146801 (2010)]. The temperature dependence of the 
Hall mobility was explained well with the conventional Bloch-Griineisen formula and yielded the 
Debye temperature 0u = 113 ± 14 K. The results indicate that the scattering of bulk electrons is 
dominated by acoustic phonons. 
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I. INTRODUCTION 


Electric and spin transport in three-dimensional (3D) 
topological insulators has received a great deal of atten¬ 
tion because of the insulators’ surface Dirac fermions 
with the unique spin texture. The presence of the 
spin texture has been confirmed by several experiments 
employing angle-resolved photoemission spectroscopy 
(ARPES), and a variety of novel transport phenomena 
have been proposed [1-8]. In real materials, however, 
unavoidable bulk conduction often hinders surface con¬ 
duction. 

Electric transport in 3D topological insulators has been 
analyzed with the so-called two-carrier model; i.e., the 
parallel conduction of carriers with different mobility [9- 
11]. Within the model, the resistivity pxx in a magnetic 
field B = pqH is expressed by the Equation 


Pxx(,B) — 


niqip,! -\- n2q2P-2 + iniqip,2 + n2q2P-i)p.iP2B'^ 


{niqipi + n2q2P2y + {niqi + 712 ^ 2 )^^?^ 2 ^^ ’ 

( 1 ) 

where qi and q 2 , ni and 122 , and pi and p 2 are respec¬ 
tively the charges, densities and mobilities of the carriers. 
Despite its usability, the model contains four active pa¬ 
rameters, which lead to large errors in most cases. The 
errors have not been discussed and their statistical re¬ 
liability remains unclear. This is a fundamental issue 
requiring a more elaborate analysis procedure. 

In this article, we report the electric transport of 
the 3D topological insulator TlBiSe 2 . The material is 
known as one of the simplest and most practical ma¬ 
terials for the transport studies [12-17]. Two-carrier 
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transport properties were precisely determined by ap¬ 
plying a newly developed analysis procedure. The pro¬ 
cedure enables an analysis with two active parameters 
and with sufficiently small errors. The magnetotransport 
properties were well explained accounting for high- and 
low-mobility electrons in the whole temperature range. 
The carrier densities were (0.61 ±0.01) x 10^® /cm^ and 
(4.45±0.01) X 10^® /cTP?, which were 10® times the values 
expected from the surface. These results indicate that 
multicarrier conduction originated from the bulk elec¬ 
trons, and the scattering was insensitive to the magnetic 
field. The temperature dependence of the Hall mobility 
exhibited a metallic behavior over the whole temperature 
range, and was explained well by the Bloch-Griineisen 
formula. This indicates that the scattering of the bulk 
electrons was dominated by acoustic phonons. It should 
be noted that analyses proposed in this article allow a 
precise identification of the two-carrier transport prop¬ 
erties and are therefore useful for all other multicarrier 
systems. 


II. EXPERIMENTAL 

The TlBiSe 2 sample used in this study was taken from 
the same batch of Ref. [16]. As reported in Ref. [16], a 
significant variation of the chemical composition was re¬ 
vealed between the starting substances and synthesized 
crystals; the variation resulted in the spontaneous elec¬ 
tron doping. The sample, which was synthesized from 
the starting composition T1 : Bi : Se = 1 : 1 : 2 in 
stoichiometric ratio, corresponded to Tli_a;Bii+a;Se 2-5 
{x = 0.064,^ = 0.29) [16, 18]. Electric transport exper¬ 
iments were performed using the conventional six-probe 
technique down to 2 K with a commercial apparatus 
(Quantum Design PPMS). Electric contacts were made 
onto the cleaved surface with room-temperature-cured 
silver paste under atmospheric conditions. Geometries 
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including the distances of electrodes and the thickness 
of the sample were determined from optical microscope 
images. Schematics of the crystal and parallel conduc¬ 
tion of the surface and bulk are presented in Figure 1(a), 
(b). The optical microscope images are also shown in 
Figure 1(c). The sample had a thickness of 0.380 mm. 
Magnetic fields (maximum strength of ±9 T) were ap¬ 
plied parallel to the c-axis. 
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FIG. 1. (Color online) (a) A schematic of the TlBiSe 2 crys¬ 
tal with hexagonal coordinates (space group: R3m [13]). 
(b) Schematics of parallel conduction of the spin-polarized 
surface and bulk. The cleaved surface forms thallium is¬ 
lands on the selenium layer [19]. (c) Optical images of the 
top surface and side surface. The sample thickness was 
d = 0.380 mm. 


III. RESULTS AND ANALYSES 


Pxx exhibited a magnetic field dependence as clearly 
seen in Figure 3(a). Here Ap^x/pxx = [Pxx{B) - 
Pxxi^)]/Pxx{B) is the magnetoresistance. For further in¬ 
sight, we transformed the two-carrier model [20]. Within 



The temperature dependence of the longitudinal resis¬ 
tivity pxx is presented in Figure 2(a). The inset shows the 
Hall resistivity pyx at 300 and 2 K. The positive tempera¬ 
ture coefficient of pxx indicates metallic conduction. The 
negative magnetic coefficients of pyx indicate electron 
conduction. These results are consistent with the ARPES 
results; the Fermi level was located in the bulk conduc¬ 
tion bands [14-16]. The pyx exhibited a linear magnetic 
field dependence at each temperature. The tempera¬ 
ture dependence of the Hall coefficient [Rh] = \Pvx/B\ 
and the Hall mobility = \Rh\/P xx are presented in 
Figure 2(b). As shown in the figure, the [Rh] exhib¬ 
ited a small temperature dependence, /th also exhibited 
a temperature variation above 15 K. [Rh] at 2 K was 
(1.775 ± 0.006) X 10“^ cm^/C and at the same tem¬ 
perature was 342±6 cm^/Vs. Because /j,h is in proportion 
to the scattering time t, the temperature variation of 
is attributed to the temperature dependence of t. We 
discuss the dependence later. 


FIG. 2. (Color online) (a) Temperature dependence of the 
longitudinal resistivity p^x in the range 2-300 K. The inset 
is the Hall resistivity pyx at 2 and 300 K. Negative magnetic 
field dependences indicate electron conduction. A linear fit of 
pyx at 2 K is also presented, (b) Temperature dependence of 
the Hall coefficient [Rh] = \pyx/PoH\ and the Hall mobility 
pn = \Rii\/pxx- |Ah| at 2 K was (1.775±0.006) x 10“^ cm^/G 
and pYi at the same temperature was 342 ± 6 cm^/Vs. 


the model, pxx is expressed by Equation (1). Equa¬ 
tion (1) was transformed to Equation (2) using n+ = 
(ni -I- n 2 ), n_ = (ni - 712 ), p+ = {pi + M2)/2, P- = 
{pi — and N = n-/rL+, M = p-/p+. In the 

present case, the charges were qi = q 2 = —e, and thus, 
—p+ > 0 and —1 < M < 1. We also defined ni > 712 ; 
i.e., 0 < Y < 1. Finally, Pxx{B), which was applied 
in further analyses, was obtained as Equation (3). The 
magnetoresistance was expressed as Equation (4). 
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Pxx{B^ - 


n+/i+ 


gi(l + N){1 + M) + 52(1 - N){1 -M) + [qi{l + N){1 - M) + 52(1 - 1V)(1 + M)](l - 


[51(1 + 7V)(1 + M) + 52(1 - N){1 - M )]2 + [ 51(1 + N) + 52(1 - 1V)]2(1 - M2)V^B2 

Pxx{i^) — 


n+e{-p+) 


1 + NM + (1 - M2)(1 - NM)^jl\B^ 


(1 + NMf + (1 - M2)2^^B2 


( 2 ) 

(3) 


{l + NMY + {l-N^M^){l-M^)p\B^ 


Within the model, Ra and pn sd H = 0, which are pre¬ 
sented in Figure 2(b), are expressed as Equations (5) 
and (6). Using /in, Equation (4) was rewritten as 
Equation (7). Because /xh was an experimental value. 


Apxx!Pxx is determined employing two additional pa¬ 
rameters, N and M. Eurthermore, pyx{B) was rewritten 
as Equation (8) using i?H and /xh, and was also deter¬ 
mined using N and M. 


i?H 


1 

n+e 


1 -h 2NM + m2 

{l + NMY 


(5) 


MH = l/i+l 


1 -t 2NM + m2 
IM/VM 


( 6 ) 


Apxx _ {l-N^)M\l-M^)plB^ _ 

Pxx (1 + 2NM + M2)2 + {1- - M2)^^B2 

[ {1 + 2NM + M^fB + {1 +NM)^{1 - 

(1 -h 2NM + M2)3 + (1 + 2NM + M2)(l - 


ApxxjPxx at 5 K and the fit with Equation 7 are pre¬ 
sented in Eigure 3(a); the fit yielded N = 0.76 ±0.01 and 
M = —0.586 ± 0.003. Apxx!Pxx at 300 K, and the curve 
expected from the N and M values with Pa at 300 K 
are also presented. The figure shows that the curve at 
300 K is consistent with the experiments. pyx{B) and 
the curves expected from the same N and M values are 
presented in Figure 3(b). The curves also agree well with 
the experimental results. We thus conclude that the N 
and M values obtained are invariant within the studied 
temperature range, and the transport properties of the 
sample are explained well with the two-channel model. 

Once the four parameters of the two-channel model 
i?H, Ph, N, and M were obtained, the densities and 
mobilities of high- and low-mobility electrons were de¬ 
termined. The densities and the mobilities at 5 K are 
summarized in Table I. Here rxi = n+(l ± N)/2, 7 x 2 = 
n+(l - N)/2, Pi = p+{l + M), and p 2 = p+{l - M). 


TABLE I. Electric transport properties at 5 K obtained from 
the two-carrier model, m and 7 x 2 , and pi and p 2 are respec¬ 
tively the densities and mobilities of the electrons. 


m 

4.45 ±0.01 

10^® /cm' 

772 

0.06 ±0.01 

10^® /cm' 

pi 

176 ±2 

cm^/Vs 

P2 

676 ±2 

cmVVs 


IV. DISCUSSION 

As demonstrated above, the magnetotransport prop¬ 
erties were well explained by the two-channel model. 
We next discuss the origin of these conduction chan¬ 
nels with the electronic band structure. A schematic 
of the structure revealed by ARPES is presented in Eig- 
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mass than that around the F point. These band char¬ 
acters are consistent with the two-carrier model, where 
the Fermi surface around the F point yields higher car¬ 
rier density and lower mobility; i.e., higher ni and lower 
/Ti. It is also concluded that the minor electrons with 
carrier density 712 , and mobility ^2 can be attributed to 
the Fermi surface around the F point. 

As discussed above, novel surface transport is hid¬ 
den in bulk metallic conduction with carrier density 
^ 10^® /eve?. In contrast, bulk insulating state with 
carrier density < 10^^ /cm^ was demonstrated for 

Tli_a;Bii+a;Se 2-5 [x < 0.028,(5 = 0.28) and the surface 
metallic conduction was distinguished [16, 17]. These 
results indicate a presence of metal-insulator transition 
in bulk between x = 0.064 and 0.028, and novel surface 
transport can be distinguished only in the bulk insulating 
samples. 
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FIG. 3. (Color online) (a) Magnetoresistance Ap^x/pxx = 
[pxx{B) — pxx{Q)]/Pxx{B) at 300 and 5 K. The fit with Equa¬ 
tion (7) at 5 K is shown by the solid line. The curve expected 
at 300 K is shown by the dashed line, (b) Hall resistivity py^ 
at 300 and 5 K. The curves expected at each temperature are 
shown by dashed lines. 


FIG. 4. (Color online) A schematic of the energy band disper¬ 
sion of TlBiSe 2 revealed by ARPES [15, 16, 21]. The smaller 
band curvature around the F point results in a larger band 
mass of electrons. The greater band curvature around the 
F point results in smaller band mass. The amplitude of the 
minimum bulk band gap is also presented. 


ure 4 [15,16, 21]. The figure shows surface Dirac fermions 
and two Fermi surfaces of the bulk; one of the Fermi 
surfaces was around the F point and the other around 
the F point. For the surface Dirac fermions, the two- 
dimensional carrier density rig = ^f/(4'^) ’"^eis estimated 
as ^ 8 X 10^^ /cm^, where fcp was the Fermi wave num¬ 
ber obtained from the ARPES spectrum (Figure 2 of 
Ref. [15]). The value of rig yielded an effective 3D carrier 
density rig/d = 2.1 x 10^^ /cm^, where d = 0.380 mm 
was the sample thickness. The estimated value of Ug/d 
was smaller than the experimental values by five orders 
of magnitude. This indicates that the contribution of the 
Dirac fermions to pxx and pyx was negligibly small; i.e., 
the transport properties are attributed to the bulk elec¬ 
trons. Figure 4 shows that the Fermi surface around the 
F point has a larger Fermi cross section and a larger band 


Finally, we discuss the temperature dependence of 
scattering realized in the system. The scattering time 
is in proportion to the mobility, and the temperature de¬ 
pendence of p,+ is thus attributed to the temperature 
dependence of average scattering time r. Note that the 
agreement with the two-carrier model indicates that the 
scattering was insensitive to the magnetic field. As dis¬ 
cussed above, the temperature dependences of N and M 
were negligibly small. These results indicate that pn can 
be expressed as pH = 0.816|/x+| using the N and M val¬ 
ues obtained (see Equation 6) and that the temperature 
dependence of /in can be attributed to r. The tempera¬ 
ture dependence of p^^ is presented in Figure 5. The fit 
with the Bloch-Griineisen formula 


1 

— oc 
r 



T 


X 


5 


dx^ 


(9) 


(e® — 1)(1 — e 
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which describes the temperature dependence of acous¬ 
tic phonon scattering, is also presented. Here 0 d is the 
Debye temperature determined in the present transport 
experiment. As shown in Figure 5, the temperature de¬ 
pendence of l//rH oc l/r agrees well with the formula, 
which yields 0 d = 113 ± 14 K. Note that 0 d agrees 
with the Debye temperature determined from the spe¬ 
cific heat [22]. The result indicates that the scattering 
was dominated by acoustic phonons. Similar results are 
often seen for metals of a weakly correlated electron sys¬ 
tem, indicating that the scattering of bulk electrons is 
explained by the conventional theory of metals. 



T{K) 


FIG. 5. (Color online) Temperature dependence of the inverse 
Hall mobility 3^1^® with the Bloch-Griineisen formula 
yields the Debye temperature 0o = 114 ± 13 K. See the text 
for details. 


V. CONCLUSION 

We investigated the transport properties of the 3D 
topological insulator TlBiSe 2 . Two-carrier transport 
properties were precisely determined by applying a newly 
developed analysis procedure. The magnetotransport 
properties were well explained accounting for high- and 
low-mobility electrons in the whole temperature range. 
This indicated that the scattering of electrons was in¬ 
sensitive to the magnetic field. The transport properties 
obtained were consistent with ARPES results, and con¬ 
tributions of the surface transport were negligibly small 
over the entire temperature range studied. The scatter¬ 
ing of bulk electrons was dominated by acoustic phonons. 
The results indicate that the bulk transport properties 
of the TlBiSe 2 sample were in the vicinity of those of 
conventional metals. These simple nature of bulk con¬ 
duction hints proper verification of novel surface trans¬ 
port properties. It should be emphasized that the analy¬ 
sis using conventional two-carrier model with four active 
parameters leads to large errors. Instead, the analysis 
proposed in this study exhibited sufficiently small errors, 
and therefore possessed statistical reliability. We also 
note that the procedure should be useful for all other 
transport studies. 
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